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Abstract: Forthcoming projects such as DES, LSST, WFIRST, and Euchd aim to 
measure weak lensing shear correlations with unprecedented precision, constraining 
the dark energy equation of state at the percent level. Reliance on photometrically- 
determined redshifts constitutes a major source of uncertainty for these surveys. 
Additionally, interpreting the weak lensing signal requires a detailed understanding 
of the nonlinear physics of gravitational collapse. We present a new analysis of the 
stringent calibration requirements for weak lensing analyses of future imaging surveys 
that addresses both photo-z uncertainty and errors in the calibration of the matter 
power spectrum. We find that when photo-z uncertainty is taken into account the 
requirements on the level of precision in the prediction for the matter power spectrum 
are more stringent than previously thought. Including degree-scale galaxy cluster- 
ing statistics in a joint analysis with weak lensing not only strengthens the survey's 
constraining power by ~ 20%, but can also have a profound impact on the calibra- 
tion demands, decreasing the degradation in dark energy constraints with matter 
power spectrum uncertainty by a factor of 2 — 5. Similarly, using galaxy clustering 
information significantly relaxes the demands on photo-z calibration. We compare 
these calibration requirements to the contemporary state-of-the-art in photometric 
redshift estimation and predictions of the power spectrum and suggest strategies to 
utilize forthcoming data optimally. 
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1. Introduction 

Weak gravitational lensing of galaxies by large-scale structure is a potentially power- 
ful cosmological probe [1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11]. Forthcoming imaging surveys 
such as the Dark Energy Survey (DES), the survey to be conducted by the Large Syn- 
optic Survey Telescope (LSST), the survey of the European Space Agency's Euclid 
satellite [12], and the proposed Wide Field Infra-Red Survey Telescope (WFIRST) 
expect to exploit measurements of weak gravitational lensing of distant source galax- 
ies to constrain the properties of the dark energy [13, 14, 15, 16, 17, 18, 19, 20, 21, 
22, 23, 24, 25, 26, 27, 28, 29, 30]. Over the last several years, it has been recognized 
that the limited precision with which the matter power spectrum can be predicted 
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may be one of several important, systematic errors that will need to be controlled in 
order to realize this goal (e.g., Refs. [31, 32, 33, 34, 35, 36, 28, 37, 38, 39]). 

A large part of this uncertainty is due to the effects of baryons on lensing power 
spectra, which were largely neglected in much of the early literature on cosmological 
weak lensing. Several groups, including our own, have begun numerical simulation 
programs designed to address this issue with large-scale numerical simulations (e.g., 
Ref. [39, 40, 41, 42, 43, 44]). With the notable exception of Ref. [33], there have not 
been detailed studies of the precision with which the matter power spectrum must be 
predicted before it becomes a relatively small contributor to the error budget. In the 
run-up to large, computationally-intensive and human resource-intensive simulation 
campaigns, we have studied the required precision with which the matter power 
spectrum must be predicted in order to realize anticipated dark energy constraints 
from cosmic weak lensing tomography. We present our results as a set of general 
guidelines on the systematic errors on the matter power spectrum as a function of 
scale and redshift. 

Another potentially dominant source of error for dark energy parameter estima- 
tors arises from the necessity of using approximate redshifts determined from photo- 
metric data rather than spectroscopic redshifts [45, 46, 47, 48, 49, 50, 51, 52, 53, 54]. 
Photometric redshifts are significantly less precise than spectroscopic redshifts, and 
can exhibit large biases. Interestingly, we find that the precision with which the 
power spectrum must be predicted is very sensitive to the precision of photometric 
redshift determinations and vice versa, a result hinted at in the Appendix of Ref. [30] . 
As part of our analysis, we model photometric redshift uncertainty and show how the 
precision with which the power spectrum must be predicted varies with photometric 
redshift errors, and conversely. 

Briefly, we find that if prior information on the photometric redshift distribu- 
tion is weak, then dark energy constraints degrade 2 — 3 times more rapidly with 
uncertainty in Ps{k) than if the photo-z distribution is characterized with high pre- 
cision. Thus we find that when photo-z uncertainty is taken into account the cal- 
ibration requirements on the theoretical prediction for the matter power spectrum 
are more stringent than previously thought. The complementarity of galaxy clus- 
tering statistics with weak lensing, well-studied in other contexts (see, for example, 
Refs. [25], [55], [56], [57], and [58]), has an ameliorating effect on power spectrum 
misestimations. Even when restricted to degree-scales, including galaxy correlation 
information can mitigate dark energy systematics induced by errors in the predic- 
tion for P5{k) by up to 50%; alternatively, neglecting galaxy clustering statistics can 
cause the statistical constraints on dark energy parameters to degrade 2 — 5 times 
more rapidly with uncertainty in either Ps{k) or the photo-z distribution. 

This manuscript is organized as follows. In §2 we describe how we model uncer- 
tainty in photometric redshifts and in Ps{k) as well as our methods for estimating 
statistical and systematic errors. We present our results in §3 and discuss their 
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implications in §4. We conclude in §5 with a summary of our primary results. 
2. Methods 

2.1 The Matter Power Spectrum 

A significant amount of the constraining power of weak lensing surveys will come from 
scales on which the nonlinear effects of gravitational collapse cannot be neglected 
[28, 59, 60]. It is possible to excise data on relatively small scales, but such an excision 
significantly degrades dark energy constraints from cosmological weak lensing [28, 59]. 
Modeling nonlinear structure formation will be essential in order for forthcoming 
galaxy imaging surveys to realize their potential for constraining dark energy and 
modified gravity [33, 61]. As this modeling is uncertain and can be an important 
source of error, we study the errors induced on dark energy parameter estimators by 
uncertainty in the underlying matter power spectrum and we quantify the relative 
importance of theoretical power spectrum errors as a function of wavenumber. Our 
results may serve as a guideline for computational programs aimed at predicting 
accurate and precise matter power spectra for the purpose of comparing with weak 
lensing data. 

In the current and past literature, the three most commonly-used techniques 
employed for predicting the matter power spectrum in the mildly nonlinear regime are 
the fitting formula of Peacock & Dodds [62], the Halo Model (see Ref. [63] for a review, 
as well as the many references therein), and the fitting formula of Smith et al. [64]. 
As a rough look at the contemporary level of uncertainty in predictions for the matter 
power spectrum Ps{k), and to set the stage for what is to come, we have plotted in 
Fig. 1 the fractional difference in Ps{k) at several different redshifts predicted by 
these three nonlinear evolution models. The differences in the predictions made by 
these methods become significant on scales {k > 0.2/iMpc~^), which coincides with 
the scales at which the constraining power of weak lensing begins to peak. When 
the baryonic physics of galaxy formation is taken into account (for example, as in 
Refs. [36, 39, 37]), among other possible effects the matter distribution within halos 
is known to change relative to N-body (dark matter-only) simulations. To illustrate 
the effect such a rearrangement may have on the matter power spectrum, in Fig. 1 
we additionally plot the fractional difference in Ps{k) that is induced when the Halo 
Model correctly predicts the nonlinear evolution but the mean concentration of dark 
matter halos (quantified by the parameter cq in Eq. 2.1) is misestimated by 20%. 

In this paper, we attempt to remain relatively agnostic about the types of errors 
that may be realized in predictions of the matter power spectrum at the percent 
and sub-percent levels, because such accurate predictions have not yet been made. 
We study two classes of models for uncertainty in the matter power spectrum. Our 
first model is motivated by simulation results suggesting that one large source of 
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Figure 1: Fractional difference in Ps{k) as predicted by the tfiree different models of 
nonlinear evolution. The superscript indicates the fiducial model, the subscript the model 
inducing the fractional change. For example, the curve labelled by represents the 
fractional change to the P<5(A;) induced by using Smith et al. to compute nonlinear power 
rather than the fiducial halo model; thus when this curve is positive the Smith et al. pre- 
diction for Ps{k) exceeds that of the halo model. The dashed curve labeled Acq shows the 
fractional change to Ps{k) induced by systematically misestimating the mean concentration 
of dark matter halos by 20%. 



systematic error in nonlinear power spectrum predictions may arise from a systematic 
error in predictions for the internal structures of dark matter halos. Significant 
rearrangement of dark matter may be a result of baryonic processes during galaxy 
formation, for example [34, 36, 37]. So, we suppose that halo abundances and halo 
clustering are well-known so that the halo model accurately predicts the gross shape 
of Pslk), but that the distribution of matter within halos is uncertain. For our 
purposes, the average halo density profile may be taken as the Navarro-Frenk- White 
(NFW) density profile [65], 

pWoc fc-^) ' (l + c-^) \ (2.1) 

V -fL200m / V -n-200m / 

The parameter c describes the concentration of the mass distribution towards the 
halo center; the mass contained within the radius -R200m defines the halo mass. The 
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average halo concentration varies with halo mass and evolves in redshift; we model 
these dependences as a power laws, 

c(m, z) = Co[m/m,,o]"(l + zf, (2.2) 

where m^fi = 2.2 x 10^^ Mq/H is the typical mass of a halo that is just starting to 
collapse at z = 0. We allow the parameters a, (3, cq to vary about their fiducial values 
of Co = 10.7, a = —0.11, and (3 = —1, so that in this model uncertainty in Ps{k) 
stems exclusively from uncertainty in the distribution of mass within gravitationally- 
coUapsed, self-bound objects. Our fiducial values for a and (3 are chosen to match 
results from the Millennium Simulation [66] , but our fiducial value of cq is roughly a 
factor of two larger than the mean concentration of Millennium Simulation redshift- 
zero halos. This larger value is motivated by the results from Ref. [36] showing that 
baryonic physics typically produces a significant enhancement to the concentrations 
of dark matter halos. We thus set our fiducial Cq according to Ref. [36] rather than 
the Millennium results. 

Our second class of model for uncertainty in the matter power spectrum Ps{k, z), 
is significantly more general. After choosing a technique for predicting a fiducial 
matter power spectrum on scales of interest, we parameterize uncertainty in Ps{k) as 
follows. The range of scales between k = 0.01/iMpc~^ and k = lO/iMpc"^ are binned 
evenly in log(/c), and uncertainty in the matter power spectrum on a scale k in the 
i^^ bin is parameterized via an additional parameter 51n(Pi), such that 

Fs{k,z)^Fsik,z)il + 6\n{P,)). 

The parameters 51n(Pi) are then allowed to freely vary about their fiducial values 
of zero. We discuss our choice for the number of S\n{Pi) parameters in §2.4. This 
scheme for studying Ps{k) uncertainty was first introduced in the context of weak 
lensing in [33]. 

2.2 Photo-z Uncertainty 

As uncertainty in photometric redshifts is likely to be one of the chief contributions 
to the error budget in future lensing measurements of dark energy, forecasts of dark 
energy parameter constraints need to include photo-z uncertainty to realistically es- 
timate the constraining power of future surveys. While it is the primary goal of this 
manuscript to study the infiuence of uncertainty in the matter power spectrum, our 
previous work has shown that photo-z calibration requirements can depend sensi- 
tively on the fiducial model for the nonlinear evolution of Ps{k)- In the Appendix of 
Ref. [30], we showed that as prior information on the photometric redshift distribu- 
tion decreases, constraints on wq degrade ~ 5 times faster if the Peacock & Dodds 
fitting formula is used to model nonlinear power rather than the fitting formula of 
Smith et al. (halo model results are intermediary between these two). The primary 
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reason for this difference stems from Smith et al. predicting more small-scale power 
relative to the other two models, a trend that is apparent in Fig. 1. These results 
demonstrate that there is a nontrivial interplay between the nonlinear evolution of 
P5{k) and the redshift distribution of the sources used to measure the weak lensing 
signal. This suggests that uncertainty in photo-z's and Ps{k) need to be treated 
simultaneously in order to accurately predict the calibration requirements for the 
matter power spectrum. 

We model the underlying redshift distribution as n{z) oc z"^ exp —{z/ zq), where 
the normalization is fixed so that n{z)dz = Na, the mean surface density of 
sources in the survey. For weak lensing studies of dark energy, the most ambitious 
planned experiments for the next ten years will be LSST and Euclid; these surveys 
correspondingly have the most stringent calibration requirements, and so we find it 
useful to quantitatively phrase our results for future very-wide-area surveys such as 
these. Thus unless stated otherwise we choose A'a = 30gal/arcmin^, and zq = 0.34, 
corresponding to a median redshift of unity. We emphasize here, though, that the 
qualitative trends in all of our results are unchanged by the particular details of the 
survey characteristics. 

We treat photometric redshift uncertainty in a relatively general manner, follow- 
ing the previous work of Ref. [45]. We assume that the source galaxies are binned 
according to photometric redshift and that the true redshift distributions of the galax- 
ies within each of the photometric redshift bins, ni{z), are related to the overall, true 
galaxy source redshift distribution, n{z), via 

miz) = niz) / " dzP^P{zP''\z), 



where zl^^ and ^^jg^ are the boundaries of the i^^ photometric redshift bin. Photo-z 
uncertainty is controlled by the function P{z^^\z). We take this to be a Gaussian at 
each redshift, 

PizP^lz) -- 



exp 



^bias) 



This may seem to be overly restrictive, but as the mean, ^bias, and dispersion, o"^, 
of this distribution can vary with redshift themselves, this parameterization allows 
for a wide variety of possible forms for the functions ni{z). We adopt a fiducial 
model for photometric redshift error in which cr^ = 0.05(1 + z), while 2bias = at 
all redshifts. We use these fiducial functions to set the values of and ^bias at 
31 control points, tabulated at intervals of Az = 0.1 between z=0 and z=3. The 
values of 2;bias and az at each of these redshifts are free parameters in our forecasts, 
so that we model photo-z uncertainty with 2 x 31 = 62 free parameters. This choice 
of binning allows for maximal degradation in dark energy constraint in the absence 
of prior information about the photometric redshift distribution of source galaxies 
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(priors that would result from, for example, a photo-z calibration program). This 
indicates that our parameterization does not enforce correlations that yield better- 
than-expected dark energy constraints. The dark energy constraints in the absence 
of prior information reduce to the same constraints that would be obtained with no 
binning in photometric redshift. 

To model uncertainty in the fiducial photo-z distribution, we introduce priors on 
the values of the dispersion and bias at the i^^ redshift control points, a* and ^^ias' 
respectively. These priors are 



where Ngp^^, is the number of spectroscopic galaxies used in each of the 31 bins of 
width Sz = 0.1 to calibrate the photo-z distribution. Our implementation of priors 
on the photo-z parameters is certainly simplistic. For example, we have further 
simplified our calculations by setting all of the N^pg^ equal to each other, effectively 
assuming that the calibrating spectra are sampled equally in redshift, whereas in 
practice there will be looser constraints on sources at high redshift than at low 
redshift. However, because the details of how a realistic calibration program will 
proceed remains uncertain at the present time, we use this simple model for prior 
information and postpone a refinement of this parameterization until the exact set 
of spectra that will be used to calibrate LSST and Euclid is better known. 

We emphasize here that Ngpec provides a convenient way to specify a one-parameter 
family of photo-z priors through Eqs. 2.4 & 2.5. The quantity Ngpec is likely not the 
true size of the spectroscopic calibration sample, but rather the equivalent size of a 
sample that fairly represents the color space distribution of the sources used in the 
lensing analysis. 

2.3 Observables 

We take the lensing power spectra of source galaxies binned by photometric redshift 
as well as the galaxy power spectra and galaxy-lensing cross spectra as observables 
that may be extracted from large-scale photometric surveys. The Limber approxima- 
tion relates the power spectrum Vs^sjik^z) associated with the correlation function 
of a pair of three-dimensional scalar fields, S; and Sj, to its two-dimensional projected 
power spectrum, Pxjxj (^) : 




(2.4) 



'bias 




(2.5) 




dz 



W,(z)W,-(;^) 



V,,,.{k = l/D^{z),z). 



(2.6) 
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Eq. 2.6 essentially describes how the two-dimensional scalar fields x; are observed as 
projections of the three-dimensional scalar fields, Sj. The angular diameter distance 
function is denoted by Da, and H{z) is the Hubble parameter. 

The weight function Wi{z) specifies the projection of the 3D source fields onto 
the 2D projected fields: 



For galaxy fiuctuations, the weight function is simply the redshift distribution of 
galaxies in the i^^ tomographic bin, ni{z), times the Hubble rate: 



The weight function associated with fiuctuations in lensing convergence is given by 



where Da(-z, z') is the angular diameter distance between z and z'. 

In principle, neither the redshift distribution of the galaxies used for the galaxy 
clustering nor the tomographic binning scheme need be the same as that used for 
cosmic shear sources, but for simplicity we use the same underlying distribution and 
binning for both so that the chief difference between the galaxy power spectrum 
Pgig^, the convergence power spectrum Pkikj, and the cross-spectrum PKigj, is the 
form of the weight functions. Above and throughout, lower-case Latin indices label 
the tomographic redshift bin of the sources. For a survey with its galaxies divided 
into Ng redshift bins used to measure the galaxy clustering, and Ng bins for the 
galaxies used to measure cosmic shear, there will be Ng{Ng + l)/2 distinct 2-D 
galaxy power spectra Pgig^, Ng^Ng + l)/2 distinct convergence power spectra Pkikj, 
and NgNg distinct cross-spectra P^gj- 

The matter power spectrum, Pg.s. (/c, z) = Vs{k, z), sources the three-dimensional 
power in cosmic shear, whereas the source power for galaxy-galaxy correlations is 
the 3-D galaxy power spectrum Vs^sj{k, z) = Vg.g.{k, z). In all of our calculations we 
restrict galaxy correlation information to low multipoles i < 300; at redshift z=l 
this corresponds to fluctuations of wavenumber k ^ 0.2/iMpc~^, so it will suffice for 
our purposes to use a linear, deterministic bias to relate the mass overdensity, S{z), 
to the galaxy overdensity, Sg{z) = b{z)6{z). 

We allow for a very general redshift-dependent bias. To model uncertainty in 
the galaxy bias function b{z), we allow the bias to vary freely about its flducial value 
of unity in A*";, galaxy bias bins, evenly spaced in true redshift, so that uncertainty 
in galaxy bias is encoded by A*";, parameters. We computed dark energy constraints 
using the Fisher analysis technique described in §2.4 for Ni, ranging from 1 to 30. 
We flnd that the dark energy constraints are insensitive to Nb ranging from 1 — 15. 




(2.7) 



Wf{z) = H{z)n;{z). 
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Throughout this manuscript, we present results pertaining to = 10 bins, so that 
the value of the galaxy bias function b{z) has independent, parametric freedom in 
redshift bins of width 6z = 0.3. While finer binning is possible, particularly if the 
number of tomographic galaxy bins Ng is increased, a further increase of parametric 
freedom is unnecessary as galaxy bias is not a rapidly varying function of redshift 
(see, for example, Ref. [67]). 

As a further simplification, we set the fiducial value of the bias function to unity 
at all redshifts, hi{z) = b{z) = 1 at all z for all i. This choice of fiducial parameter 
values is conservative, because the galaxies observed as part of high redshift samples 
will likely be biased [68], exhibiting relatively stronger correlations than matter, so 
we underestimate signal-to-noise of galaxy clustering measurements in the fiducial 
case. 

2.4 Parameter Forecasting 

We estimate the constraints from upcoming photometric surveys using the formalism 
of the Fisher information matrix. Useful references for this formalism include [69, 
70, 71, 72, 24]. The Fisher matrix is defined as 



The parameters of the model are pa and the Oa are the observables described in 
§ 2.3. Greek indices label model parameters while Latin, upper-case indices label 
distinct observables. We take imm = 2 for all observables. For the lensing spectra, 
we take P^k observables we set f^ax = 3000 so that the assumptions of weak lensing 
and Gaussian statistics remain relatively reliable [73, 74, 75, 76, 6]. For the galaxy 
clustering statistics, we eliminate small-scale information so that we do not need 
to model scale-dependent galaxy bias, which is potentially complicated in itself. 
Therefore, we set £max = 300 for Pgg and P^g, corresponding roughly to angular 
scales of ~ 1 degree. We emphasize that this restriction is very conservative as it 
implies that our joint analysis does not employ the use of Baryon Acoustic Oscillation 
features in our galaxy power spectra. As we will see, the added benefit of a joint 
analysis stems primarily from the increased ability to self-calibrate parameterized 
uncertainty in Ps{k) and the distribution of photometric redshifts. 

In Eq. 2.8, C^g is the inverse of the covariance matrix; our treatment of the 
covariance matrix calculation, and its associated Fisher matrix, is very similar to 
that in Ref. [77], to which we refer the reader for additional details. Briefly, the 
covariance between a pair of power spectra P^.^. and -Px^xn is given by 






(2.9) 
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where in the case of either galaxy power or convergence power the observed spectra 
PxiXj have a contribution from both signal and shot noise, 

where N^.g. = ^yN^ is the shot noise term for galaxy spectra, with Nj^ denoting the 
surface density of sources, and N^iKj = 5ij7intNf' is the shot noise for convergence. 
We calculate the observed cross-spectra P^igj without a contribution from shot noise, 
so that Pi^.g. = -Psigj, because galaxies are lensed by mass separated from them by 
cosmological distances, so the cross-correlation of the noise terms should be small. We 
adopt a common convention of setting the intrinsic galaxy shape noise 7int = 0.2 and 
absorb differences in shape noise between different surveys into the surface density 
of sources A^a- 

The inverse of the Fisher matrix is an estimate of the parameter covariance near 
the maximum of the likelihood, i.e. at the fiducial values of the parameters. There- 
fore, the measurement error on parameter a marginalized over all other parameters 
is 

aipa) = V[F-%a. (2.10) 

Gaussian priors on the parameters are incorporated into the Fisher analysis via F^^ 
in Eq. 2.8. If one is instead interested in unmarginalized errors, for example to test 
the intrinsic sensitivity of the observables to a parameter pa, the quantity y^l/[F]aa 
provides an estimate of the uncertainty on pa in the limit of zero covariance between 
Pa and any of the other parameters in the analysis. 

The Fisher formalism can also be used to estimate the magnitude of a bias that 
would occur in parameter inference due to a systematic error in the observables. 
If AO A denotes the difference between the fiducial observables and the observables 
perturbed by the presence of the systematic error, then the systematic offset in the 
inferred value of the parameters caused by the error can be estimated as 

= J2[F-%,J2^2i + l)U,Y.^^^^AB^- (2-11) 

13 I A,B 

We assume a standard, flat ACDM cosmological model and vary seven cosmolog- 
ical parameters with fiducial values are as follows: VLm^? = 0.13, wq = —1, = 0, 
Qbh^ = 0.0223, ns = 0.96, ln{AD = -19.953, and = 0.73. We utilize the following 
marginalized priors: A^m^^ = 0.007, AVl^h^ = 0.001, An^ = 0.04, Aln(A|) = 0.1. 
These priors are comparable to contemporary uncertainty [78], so this choice should 
be conservative. We have verified that strengthening these priors to levels of uncer- 
tainty that will be provided by Planck [79] does not induce a significant change to 
any of our results. 

We determined the number of independent parameters for the matter power 
spectrum with an analysis of the off-diagonal elements of the inverse Fisher matrix. 
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The parameter covariance is 




■-1 



(2.12) 



In general, —1 < Qais < 1, with Qa^ = (— )1 corresponding to the case where 
parameters Pa and are perfectly (anti-)correlated and Qa^ = corresponding 
to uncorrelated parameters. By increasing the number of matter power spectrum 
parameters until the Fisher matrix is no longer invertible in the absence of prior 
information on these parameters, we determined that a lensing-only analysis reaches 
a level of total information loss when ten 6\n{P[) parameters are used. By studying 
the behavior of the off-diagonal Fisher Matrix entries as the number of matter power 
spectrum parameters are increased, we find that this state of information loss occurs 
after values of Qa/s = ±0.8 obtain between pairs of distinct 51n(Pi), in agreement 
with the method used in Ref. [33] to arrive at this conclusion (D. Huterer, private 
communication). Because galaxy correlation observables provide additional informa- 
tion with which to self-calibrate matter power spectrum parameters, including Pgg 
and P^g allows for slightly finer binning in wavenumber, but for the sake of facili- 
tating a direct comparison between the different sets of observables we have limited 
our analysis to ten parameters (51n(Pi), irrespective of whether we consider a joint 
analysis or weak lensing alone. 

There is additional freedom in the choice of the number of tomographic bins one 
uses to divide both the sources used to measure lensing as well as the sources used 
to measure galaxy correlations. As has been noted in previous studies, for example 
Ref. [45] , dark energy information from lensing saturates at iV^ = 5 tomographic bins; 
this saturation point is determined by using the Fisher matrix (in the limit of perfect 
prior knowledge of all nuisance parameters, in our case the photo-z parameters o"* and 
2;^ias! and the (51n(Pi) parameters) to compute the statistical constraints a{wo, Wg.) and 
increasing the number of tomographic bins until the constraints cease to improve. 
We find that for the case of galaxy clustering this information saturation occurs at 
Ng = 10 tomographic bins, although we note that this saturation point depends on 
the maximum multipole used in the analysis. For example, Ref. [55] uses i^ax = 2000 
for their galaxy clustering analysis and finds that additional information is available 
by increasing the number of bins to A^^^ = 40. Our smaller Ng saturation point is a 
consequence of our conservative choice for i^ax, and the lack of BAO information 
implied by this choice. 

3. Results 

3.1 Power Spectrum Self-Calibration 

We begin by presenting our calculation of the scale-dependence of the sensitivity of 
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weak lensing (with and without galaxy correlations) to the matter power spectrum. 
Our second model for Ps{k) uncertainty is well-suited to this investigation: the 
constraints on parameter 51n(Pi) provide an estimate of the statistical significance of 
the weak lensing signal produced by correlations in the matter distribution on scales 
k ^ ki. In the prevaihng jargon, this calculation corresponds to self calibration of the 
matter power spectrum. 

The constraints from this computation are plotted as a function of scale in Fig. 2. 
The magenta curves at the bottom of the Fig. 2 pertain to unmarginalized constraints 
on the (51n(Pi). In other words, covariant uncertainty in cosmology, photometric red- 
shifts, and galaxy bias is not taken into account in the magenta curves. In plotting 
the red and blue curves we illustrate our results when this covariance is accounted 
for by marginalizing over all other parameters in the analysis. The red curves corre- 
spond to a calculation with Ngpec = 8000, or Acr^/a^ ~ 10~^. The blue curves pertain 
to a calculation with Ngpcc = 2 x 10^, which is sufficiently large that further increases 
to Ngpec do not improve constraints on any of the parameters in our analysis, so 
priors this tight effectively correspond to the case where the photo-z distribution is 
known perfectly. The solid curves include galaxy correlation observables (Pgg and 
Pi^g) in addition to lensing observables and thus lie strictly below the dashed curves 
(lensing only). The step-like appearance of the curves reflects the coarse binning in 
wavenumber of our parameterization of Ps{k) uncertainty: a forthcoming very- wide- 
area, LSST- or Euclid-like survey is only able to constrain ~ 10 independent matter 
power spectrum parameters (see § 2.4). While each of the curves in Fig. 2 pertains 
to a calculation in which we used the Smith et al. [64] fitting formula as our fiducial 
Ps{k)., the results using either the halo model or the Peacock & Dodds [62] fitting 
formula are nearly identical, so conclusions drawn from Fig. 2 are quite robust to 
detailed changes in the fiducial model for nonlinear collapse. 

The minimum of the unmarginalized constraints in Fig. 2 at /c ~ 2/iMpc~^ 
occurs on the scale at which weak lensing is most intrinsically sensitive to matter 
overdensities. This minimum occurs on a physical scale nearly an order of magnitude 
smaller than the minimum of the marginalized constraints. This observation is an 
extension of the previous work of Ref. [33] and is itself an important result as it 
demonstrates the need for precision in the prediction for Ps{k) over the full range of 
nonlinear scales k < 5 /iMpc~^. Because this shift in scale-dependence occurs even 
for the case of perfect prior knowledge on the photo-z parameters (Ngpcc = 2 x 10''), 
then it is not the effect of photometric redshift uncertainty that drives this shift in 
scale, but rather degeneracy with cosmological parameters. An analysis of the off- 
diagonal Fisher matrix elements shows that covariance of the 6ln{Pi) with dark energy 
parameters is chiefly responsible for this dramatic shift in the scale-dependence. To 
be speciflc, with = (51n(Pi) and pj = WQ,Wa., the corresponding (Eq. 2.12) have 
the maximum magnitudes of any of the cosmological parameters in our parameter 
set and they attain their maxima at ~ 0.1/iMpc~^. 
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Figure 2: Statistical constraints on 5ln(P{), the parameters encoding uncertainty in the 
cahbration of the matter power spectrum, are plotted against the scale of the wavenumber. 
We plot unmarginalized constraints as magenta curves; the minimum of the magenta curves 
at « 2 /iMpc~^ illustrates that the intrinsic sensitivity of the weak lensing signal peaks 
at this scale. With the red and blue curves we plot marginalized constraints on the (51n(Pi) 
parameters for different levels of uncertainty in the distribution of photometric redshifts. 
The vertical axis values for the red and blue curves give the statistical precision with which 
a future very-wide-area survey such as LSST or Euclid will be able to self-calibrate the 
theoretical prediction for the matter power spectrum on a given scale. 



3.2 Statistical constraints on Dark Energy 

We proceed with results on the sensitivity of dark energy constraints to uncertainty 
in predictions of the nonlinear evolution of Ps{k), incorporating possible additional 
uncertainty from photometric redshift errors. Results for the Sln{Pi) model appear in 
3.2.1 while those pertaining to the more restrictive 'halo model" treatment of power 
spectrum uncertainty are discussed in 3.2.2. 

3.2.1 The (51n(Pi) Model 

In Figure 3 we depict contours of the degradation in the statistical constraints on 
the dark energy equation of state associated with simultaneous uncertainty in Ps{k) 
and photometric redshifts. The constraints on wq and Wa are shown in units of the 
perfectly calibrated limit, when power spectra and photometric reshifts are known 
so well as to be inconsequential to the dark energy error budget; we denote the 
constraints on wq and in the limit of perfect calibration by a^^^^{wQ) and a^^^^{wg), 
respectively. For example, we plot the ratio S = cr{wo)/(r^^^^{wo) to illustrate the 
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Figure 3: Contours of degradation in wq and constraints appear in the top and bottom 
panels, respectively. The degradation is quantified by H = a{wo,Wa)/<7^'^'^^{wo,Wa.), where 
cjP*^''^ is the level of statistical uncertainty of an LSST- or Euclid-like survey in the limit 
of perfect certainty on the photo-z distribution and the nonlinear evolution of Ps{k). The 
precision in the calibration of the matter power spectrum Ps{k) appears on the vertical 
axes, while the priors on the photo-z distribution, parameterized by Ngpec, appear on the 
horizontal axes. Prior information about the functions Zbias and az, which govern the 
uncertainty in the photo-z distribution P{zP^\z), is distributed uniformly in redshift ac- 
cording to Azbias = \/2A(T^ = az/ Y^Nspec- Priors on the 6ln{Pi) parameters are distributed 
uniformly in among bandpowers in log(fc). Dashed curves pertain to a survey using only 
weak lensing information, solid curves to a joint analysis that includes galaxy clustering. 
The gray, horizontal lines roughly bound the range of matter power spectrum uncertainty 
that is attainable in advance of these surveys. 
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Table 1: Baseline Constraints 



Observables 




Weak Lensing Only 0.071 
Joint Analysis 0.058 



0.22 
0.19 



0.022 
0.018 



Notes. — Column (1) specifies whether or not observables employing galaxy clus- 
tering (Pgg and P^g) were used in the calculation. Columns (2), (3), and (4) give 
the statistical constraints on wq, w^, and Wp^, respectively, that can be obtained by 
future very- wide-area surveys such as LSST or Euclid. Note that these constraints 
account for statistical errors only and so represent the optimistic limit of achievable 
dark energy constraints for a survey with these characteristics. 

level of degradation of wq constraints. For the sake of scaling our results to an 
absolute statistical constraint, we summarize these baseline constraints in Table 1. 
We refer the reader to § 5 for a discussion of how the results we present here change 
when considering a survey with characteristics similar to DES. 

The levels of constraint degradation depend upon the precision of both the matter 
power spectrum, parameterized by 51n(Pi), and the photometric redshift distributions 
of sources. We quantify the precision of power spectrum prediction by a prior con- 
straint on the 51n(Pi) parameters, A6ln{P]). For example, a value of A51n(Pi) = 0.1 
corresponds to a 10% precision on the bandpower at wavenumber ki. For simplicity, 
we apply the same prior at all wave bands in order to produce Fig. 3. The assumed 
level of calibration of photometric redshifts is specified by the parameter Nspec, as 
discussed in § 2.2. 

The horizontal gray lines roughly bound the range of precision in the prediction 
for Ps{k) on scales relevant to lensing and large-scale galaxy clustering that may 
be attainable by near-future numerical simulation campaigns. The precise values 
of the realized precisions will depend upon the resources dedicated to address this 
issue as well as the ability of data to constrain baryonic processes that alter power 
spectra and numerical simulations to treat these baryonic processes. For example, 
the Coyote Universe simulation campaign [42] has already achieved a 1% calibration 
of the matter power spectrum to scales as small as /c ~ 1 hMpc~^] future results 
from, e.g., the Roadrunner Universe [80] will improve upon these results, although 
it is still not clear in detail how precisely Ps{k) will be calibrated to scales as small 
as /c ~ 10 /iMpc~^, especially when uncertainty in baryonic physics is taken into 
account. 

Fig. 3 contains contours for parameter degradation when using lensing data alone 
(dashed curves) as well as the corresponding constraint degradations when both 
lensing and galaxy clustering observables are considered (solid curves). To be sure. 
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the additional information available to a joint analysis guarantees that a^'^^^{wQ,Wa) 
from the joint galaxy clustering and lensing analysis is less than the corresponding 
constraint when considering lensing observables alone. In each case, we plot the 
degradation S = a{wo,Wa)/o'^'^^^{wo,Wa) relative to the idealized constraints for that 
technique. Accounting for this difference, Fig. 3 highlights the dramatic relaxation 
on the calibration requirements of both power spectrum predictions and photomet- 
ric redshifts provided by including galaxy clustering statistics. For the purpose of 
studying dark energy, this relaxation is the most significant advantage provided by 
utilizing galaxy correlations in the analysis. 

To further illustrate this point, consider particular examples that can be gleaned 
from Fig. 3. In the limit of perfect prior knowledge of Ps{k) and the photo-z distri- 
bution, the increase in constraining power that a joint analysis has over an analysis 
that includes lensing observables only is merely ~ 15 — 20%; however, as the preci- 
sion in the calibration of both photometric redshifts and the matter power spectrum 
decreases, the dark energy constraints degrade by a factor of 2 — 5 more rapidly 
when galaxy clustering information is neglected. For example, suppose that pre- 
dictions for the matter power spectrum attain 1% level of precision in advance of 
LSST or Euclid. A joint analysis with the statistical equivalent of Nspoc ~ 20, 000 (a 
realistic proposition as we are assuming these calibrating galaxies to be distributed 
evenly in redshift) yields constraints on Wg, that are weakened by 75% due to pho- 
tometric redshift uncertainty, whereas if galaxy clustering information is neglected 
Nspcc ~ 100,000 will be required to protect against the same level of degradation in 
dark energy constraints. 

3.2.2 The Halo Model 

In Figure 4 we display the degradation in the statistical constraints on wq (left panels) 
and Wa (right panels) due to uncertainty in the halo concentration parameters cq, 
a, and /3 [Eq. (2.2)]. The precision level of the calibration of cq appears as the 
horizontal axis in the top panels, a in the middle panels, and /3 the bottom panels. 
Along each row of panels in Fig. 4, only the parameter labeled on the horizontal axes 
has uncertain prior knowledge. In other words, the remaining two concentration 
parameters are treated as perfectly known for simplicity. As in Fig. 3, each curve is 
normalized to the constraint that would be realized in the limit of perfect knowledge 
of the halo structure parameters, a'^^'^'^^. Note however, that this baseline constraint 
has been recomputed for each assumed value of photometric redshift uncertainty, 
as specified by the color-coding of each curve, so that the degradation represents 
only that amount of additional degradation due to uncertainty about halo structure. 
Thus for each level of photometric calibration precision appearing in the legend, 
Nspec = 5000, 10^, 10^, and 10®, the vertical axis value gives the degradation in 
the Wo{ws) constraints strictly due to uncertainty in the parameter labeled on the 
horizontal axis. 
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Figure 4: Plot of the degradation in dark energy constraints as a function of priors on halo 
concentration parameters, defined by c(m, z) = co[m/m*^o]"(l + z)^. Degradation of the 
constraints on wq are shown in the left panels, while degradation of constraints are shown 
in the right panels. The top panels show degradation as a function of the fractional prior 
uncertainty in the normalization of the mass concentration relation cq. The middle row of 
panels shows degradation as a function of the power law index describing the dependence 
of concentration on mass, a, and the bottom panels show the degradation as a function of 
the power-law index describing the dependence of concentration on redshift, f3. Different 
levels of priors on photo-z parameters are color-coded as labeled in the middle panels. 



In all cases, the degradation of dark energy parameters is relatively modest. In 
particular, the degradation induced by halo structure uncertainty alone is < 15% 
for all reasonable models, and significantly less if photometric redshifts are well cali- 
brated or if galaxy clustering statistics are employed. This is in qualitative agreement 
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with Ref. [28], who studied cahbrating halo structure parameters as a means to ac- 
count for the influence of baryonic processes on the lensing power spectrum. It is 
not expected that halo structure alone can account for all of the effects of baryonic 
processes on the power spectrum (e.g., Refs. [32, 36, 28, 39]). However, combining 
the results of Fig. 3 with Fig. 4 suggests that if uncertainty in the power spectrum 
due to baryonic processes can be modeled by concentrations with a residual of order 
~ 1%, then dark energy parameter degradation induced by such uncertainty may be 
limited to quite modest values. It may even be possible to incorporate additional 
parameters to describe, for example, the hot gas components of groups and clusters 
[32, 36, 39], at a relatively modest statistical cost. 

3.3 Systematic Errors on Dark Energy 

In this section we explore the related problem of systematic errors on the dark energy 
equation of state parameters induced by uncertainty in predictions of the matter 
power spectrum. To this point, we have already addressed degradation in dark 
energy parameters induced by treating uncertainty in the matter power spectrum as 
a statistical uncertainty. In such a calculation, the underlying assumption is that a 
model for the power spectrum is accurate but the parameters of the model are known 
with imperfect precision. Our aim in this section is to elucidate dark energy equation 
of state errors in the related circumstance of a systematic error on power spectrum 
predictions. The underlying framework is that a model for the power spectrum exists 
and is assumed to be correct (or at least, that it contains the true power spectrum 
within its parameter set), but the model misestimates the power spectrum over some 
range of wavenumbers. This situation may be the most relevant to forthcoming data 
analyses if, as a specific example, the dominant errors in power spectrum predictions 
stem from systematic errors in the numerical treatment of baryons. The result of 
parameter inferences that use such systematically-offset theoretical power spectra will 
be systematically-offset dark energy equation of state estimators. In this section, we 
quantify the systematic errors in dark energy equation of state parameters induced 
by systematic offsets in the matter power spectrum. Not surprisingly, the results of 
the previous section will be a useful guide for anticipating and understanding the 
results of this section. As before, we begin with our more general 6ln[Pi) model and 
later address systematic errors in the halo-based model. 

3.3.1 The (51n(Pi) Model 

We begin our treatment of systematic errors with our 6ln{P[) model. This model 
treats the predicted power spectrum as a sequence of bandpowers distributed evenly 
in log(/c). The generality of this model enables us to address a specific, important 
concern, namely: How does the severity of the induced systematic error in the dark 
energy equation of state depend on the comoving scale k at which the prediction for 
Ps{k) is erroneous? In so doing, we can prescribe how well the theoretical matter 
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power spectrum should be cahbrated as a function of scale, setting a specific goal for 
large-scale simulation efforts. 

In Fig. 5 we have plotted the systematic error induced on wq (top panel) and 
Wa (bottom panel) by a 5% error in Ps{k) at the comoving scale k labeled along the 
horizontal axis. The induced systematic errors have been calculated according to 
Eq. 2.11, where the systematic shift to the observables is induced by introducing a 
5% systematic error in Ps{k) over a range of wavenumbers with a width that is 10% 
of ken, the scale at which the error is centered, viz. 



The curves in Fig. 5 are color-coded according to the priors on the photo-z dis- 
tribution, and the magnitude of the induced systematic error has been normalized 
by a'^P'^''^, the statistical uncertainty of the parameter assuming perfect prior knowl- 
edge of each of the S\n{Pi) but with photo-z uncertainty at the level given by the 
color-coding of the curve and the observables included in the analysis. 

Each curve approaches zero at large and small scales. Very large scale lensing 
correlations contribute little constraining power on dark energy, and so it should be 
expected that scales larger than k < 10~^ /iMpc"^ should contribute comparably 
little to the error budget; errors in Psik) on scales smaller than k > 5/iMpc~^ are 
comparably tolerable because we make no use of correlation information from multi- 
poles greater than i^ax = 3000. From the wavenumber-range containing the peaks in 
curves in Fig. 5 we can see that the most significant dark energy biases induced by 
systematic errors in Ps{k) come from errors on scales 0.05 /iMpc"^ ^ ^ ^ 3 hMpc"^. 
This finding is consistent with the results presented in Fig. 2, in which we can see 
that this is the same range of scales over which the tightest statistical constraints 
can be obtained by self-calibrating the parameters Sln{Pi). Both of these figures thus 
illustrate that the weak lensing information about dark energy that will be available 
to future very-wide-area surveys such as LSST or Euclid comes primarily from grav- 
itational lensing events produced by perturbations on scales 0.05 /iMpc"^ ^ ^ ^ 3 
/iMpc"^, in good agreement with previous results [33]. 

Although almost all of the curves are limited to systematic biases (5(wo,Wa) ^ 
l(T(wo,Wa), we remind the reader that the y-axis value gives the magnitude of the 
bias when the systematic error in Ps{k) is isolated to just a single bin of wavenumbers 
of width Ak/k = 0.1, and when the magnitude of the power spectrum error is 5%. 
However, we only chose these particular values for the sake of making a definite 
illustration, and so for the results in Fig. 5 to be useful to the calibration program 
it will be necessary to scale the systematics we predict according to the particular 
details of the matter power spectrum error whose consequences are being estimated. 
We give several examples of this below to illustrate the utility of our calculations. 




(3.1) 
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Figure 5: Plot of the systematic error in the inferred value of wq (top panel) and Wg, 
(bottom panel) induced by a systematic misestimation of the matter power spectrum Ps{k) 
over a small range of wavenumbers. The value of the induced systematic is scaled in units 
of the statistical uncertainty of an LSST- or Euclid-like survey, and is plotted against the 
scale at which Psik) has been incorrectly predicted by 5% over a range of wavenumbers 
of width Ak/k = 0.1. See the text for a detailed description of how these results can be 
rescaled for Ps{k) errors spanning a range of wavenumbers of a different width. Results 
pertaining to an analysis that only uses weak lensing information appear as dashed curves, 
a joint analysis that includes galaxy clustering as solid curves. We show results for two 
different levels of photo-z calibration, with the red curves corresponding to Ngpec = 10'^, 



and the blue curves to N, 



spec 



Suppose there is a 3% error in Ps{k) made over a range of Ak/k = 0.4, centered 
at A; ~ 2 hMpc~^. For definiteness, consider a lensing-only analysis with photo-z 
uncertainty modeled by Ngpec = 10^. The y-axis value of the corresponding (dashed, 
red) curve at A; ~ 2 hMpc~^ is Swq = 0.9a{wo)- This value needs to be rescaled by 
a value of 3/5 to account for the difference between the magnitude of this example's 
actual error in Ps{k) and the 5% error plotted in Fig. 5; additionally, scaling by a 
factor of 4 is necessary to account for the fact that the range of scales over which 
the error is operative spans 4 of our bins^, giving an estimate of 6wo ~ 2.2a{wo)- 

If a 3% Ps{k) error, again spanning 4 of our Ak/k = 0.1 bins, is instead made 
at A; ~ 0.8 /iMpc"^, the estimation method of the first example naively implies that 
there is zero systematic error associated with such a power spectrum misestimation 
because the error is centered at a wavenumber at which Swq changes sign, and so 
the contributions to the net dark energy systematic to the left and right of /c ~ 0.8 
/iMpc~^ appear to cancel. Such cancellations are not necessarily spurious, and in 
fact a very general formalism for choosing a set of nuisance parameters specifically 
designed to take advantage of this phenomenon has recently been proposed [81]. 
However, because it may not be known if the sign of the matter power spectrum 
error inducing the biases also changes sign over the range of wavenumbers on which 
the error is made, it may not be possible to exploit this sign change to minimize the 
net effect of the error. In such a case, we advocate assuming the worst-case scenario, 
that biases produced by errors in multiple bins in wavenumber conspire to contribute 
additively; by construction this will yield a conservative estimate for the systematic 
induced by the error in Ps{k). Thus for this example, the absolute value of the y-axis 
values should be used to estimate the net dark energy systematic. The magnitude of 
the dashed, red curve for this particular systematic peaks at Swq = 0.5cr(wo) at the 
endpoints of its operative range, 0.6 hMpc-^ <k<l /iMpc"^; the curve approaches 
these maxima from its value of zero at A; = 0.8 /iMpc"^ and so we approximate this 
as an 0.25a{wo) error spanning 4 of our bins. Thus for a 3% systematic error in 
Psik) spanning the range 0.6 /iMpc~^ ^ k ^ 1 /iMpc"^, our final estimate is given 
by Swo = [0.25 X 4 X (3/5)] a{wo) = 0.6a(u;o)- 

For power spectrum errors made over very broad ranges of wavenumber, the 
simple linear scaling of Ak/k is no longer appropriate and one must rely on the full 
machinery of our calculation to integrate the absolute value of 6wo/(t{wo) over the 
scales over which the error is operative. This is also useful to conservatively estimate 
an ultimate target goal for the matter power spectrum calibration effort. When per- 

^ Simple, linear scaling is a very good approximation when correcting for the magnitude of the 
Psik) error, but this prescription is only approximately correct when rescaling according to the 
width of the range of scales over which the error is made. We find that the derivatives of lensing 
power spectra with respect to 6\n{Pi) parameters are stable to roughly factor-of-five changes in 
numerical step-size, and so simple, linear scaling will be appropriate so long as the width of the 
wavenumber range is less than Ak/k < 0.5. 



- 21 - 



forming this integration on scales l/iMpc^ ^ k < S/iMpc^ , the worst-case estimate 
is a systematic error of 5(wo) ~ 3 — 4(t(wo);^ if one instead integrates a 5% Ps{k) error 
over the entire range of wavenumbers 0.01/iMpc~^ ^ k ^ 5/iMpc~^, a worst-case es- 
timate of the coherently contributing biases ranges from 9 — 13cr. To ensure that dark 
energy systematics are kept at or below the level of statistical constraints, 0.5% accu- 
racy in the prediction for Ps{k) over the entire range of 0.01/iMpc"^ < A; < 5/iMpc"^ 
will be required of the simulations calibrating the matter power spectrum. Note that 
these requirements are somewhat more restrictive than estimations from previous 
work [33]. 

3.3.2 Dependence on Multipole Range 

The results presented above in §3.3.1 depend sensitively on imax, the maximum mul- 
tipole used in the cosmic shear analysis. Naturally, as £max increases the matter 
power spectrum must be modeled with greater precision and to smaller scales. Fig- 
ure 6 represents a simple illustration of this point. Each curve in Fig. 6 pertains 
to a cosmic shear-only experiment in the limit of perfect knowledge of the photo-z 
distribution, but with different choices for imax color-coded according to the legend. 
The axes are the same as those in the top panel of Fig. 5, and so this plot shows how 
wq biases induced by Ps{k) errors change with the choice for the maximum multipole 
used in the lensing analysis. 

At the end of §3.3.1 we described how to use results such as those appearing 
in Fig. 6 to estimate the precision to which the matter power spectrum must be 
predicted in order to guarantee that dark energy biases induced by Ps{k) systematics 
are kept at or below the level of the statistical constraints. Briefly, one assumes the 
worst case scenario, that the sign of the Ps{k) errors conspire to contribute coherently 
to the dark energy bias. In this case, to estimate the most severe dark energy bias 
that could be induced by an error in Ps{k) made over a range of wavenumbers, 
one simply adds the absolute value of the relevant curve over the relevant range of 
wavenumbers. To obtain a more optimistic estimation, one could suppose that the 
errors in each bin are perfectly uncorrelated, in which case they may be treated as 
independent Gaussian random variables so that their net contribution to the error 
budget is computed by adding the individual contributions in quadrature. In either 
the pessimistic or optimistic case, one obtains a precision requirement by finding the 
magnitude of the Pslk) error that would result in the sum described above equal 
to unity, since this would imply that the net systematic bias on the dark energy 
parameter is equal to the statistical constraint on that parameter (recall that the 
curves in Figures 5 and 6 plot the systematics in units of the statistical uncertainty 
of the survey). Ideally, of course, the goal of the Ps{k) calibration program is to 
achieve sufficient precision such that the systematics are well below the level of 

^The exact number depends on the level of photo-z calibration as well as the choice of observables. 
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Figure 6: Plot of the systematic error in the inferred value of wq (top panel) induced by 
a systematic misestimation of the matter power spectrum Fsik). This plot is identical to 
Figure 5 except the curves illustrate weak lensing-only results in the limit of perfect core 
calibration (Ngpec — ?• oo) but with different choices for the maximum multipole. Naturally, 
the sensitivity of the experiment to small-scale Psik) errors decreases as ^max decreases. 



statistical constraints, and so setting this sum to unity simply provides a guideline 
for the calibration. 

In Fig. 7 we have performed the calculation of the power spectrum precision 
requirement as a function of ^max- AH curves pertain to weak lensing-only experi- 
ments^, with the level of photo-z precision coded according to the legend. The top 
three curves correspond to the optimistic precision estimate (errors in each bin are 
added in quadrature), the bottom three curves to the pessimistic estimate (the abso- 
lute value of the error in each bin are added). The requirements plotted in Fig. 7 are 
set according to the magnitude of the systematics in wq; the Wa-based requirements 
are very similar. 

The optimistic calculation assumes that the Ps{k) errors at different wavenumber 
are completely independent; the pessimistic calculation assumes that the errors are 
perfectly correlated. Because the level of correlation between matter power spectrum 
errors at different wavenumbers will not be known, we stress that for a given i^ax 
the only way to guarantee that Ps{k) errors do not contribute significantly to the 
dark energy error budget is to attain the level of precision illustrated by the bottom 
curves. 

^We have not illustrated the significance of including galaxy clustering statistics because it is a 
relatively minor effect, as evidenced by Fig. 5. 
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Figure 7: Plot of the required precision with which the matter power spectrum must 
be predicted as a function of £max- The top curves represent an optimistic estimate, the 
bottom curves a pessimistic estimate. Ah curves represent results for a cosmic shear-only 
analysis; calculations for different levels of photo-z precision are color-coded according to 
the legend. To guarantee that matter power spectrum errors do not contribute significantly 
to the dark energy error budget, the precision in the prediction for Psi^) niust reach the 
level illustrated by the bottom curves. One uses the results in Figure 8 to estimate the 
wavenumber to which the level of precision plotted here must be attained. See text for 
details concerning the calculation of these estimates. 



In Fig. 7, the crossing of the Ps{k) precision requirement curves pertaining to 
different levels of Ngpec may seem somewhat counterintuitive; one might expect that 
improving photo-z uncertainty can only lead to more stringent demands on the ac- 
curacy of the prediction for the matter power spectrum. This intuitive expectation 
is supported by a theorem proved Appendix A of Ref. [53], in which the authors 
demonstrate that the net Ax^ induced by a systematic error is always reduced by 
the addition of (unbiased) prior information. However, as shown in Appendix B of 
the same paper when marginalizing over multiple parameters, Ax^ per degree of free- 
dom may increase. In our case, for large maximum multipoles (^max ^ 1000) adding 
prior photo-z information and marginalizing over our 62 photo-z parameters leads to 
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Figure 8: Plot of k^^, the wavenumber to which the matter power spectrum must be 
predicted, as a function of ^max- The quantity fe^q^ is insensitive to the level of photo-z 
precision, and so the same pertains to all values of Ngpec- One uses the results of 

Figure 7 to estimate the level of precision in the prediction for Psik) that needs to be 
attained for all k < k^ff^. 



a mild increase in the Wq bias, a fact which we have traced to a mild difference in 
the degeneracy between wq and for different values of Ngpec- 

From Fig. 6 it is evident that for smaller choices of £max dark energy biases are 
less sensitive to matter power spectrum errors on small scales. For each i^nax there 
is a maximum wavenumber, fc™q^, such that systematic errors in Psik) for k > k'^^ 
do not produce significant biases in dark energy parameters. We estimate k^^ by 
finding the bin in wavenumber at which the induced wq bias becomes less than 10% 
of the maximum magnitude that 6wo attains in any bin. Because of the steepness of 
the scaling of 6{wo)/a{wQ) with k on scales smaller than the wavenumber at which 
the systematics attain their maximum magnitude, we find that our k^^ estimations 
are insensitive to the choice for this percentage. We present our results for k^^ as a 
function of £max in Fig. 8. 

The results in this section provide a set of concrete benchmarks for the campaign 
of numerical simulations designed to calibrate the prediction for the matter power 
spectrum, as well as a guideline for choosing the maximum multipole that should 



be included in any cosmic shear analysis. For a given ^max, one uses the results 
presented in Fig. 7 to estimate the precision with which Ps{k) must be predicted on 
all scales k < k^^, where the k^^ estimate appears in Fig. 8. Of course choosing 
smaller values of £max naturally decreases the constraining power of the survey, and 
so the results we present here can be used to inform the optimal choice for imax 
that balances the need for statistical precision against the threat of matter power 
spectrum systematics. 

3.3.3 Halo Model 

In Fig. 9 we illustrate our results for the propagation of systematic errors in halo 
concentration parameters through to Wq (top panel) and (bottom panel). Each 
curve corresponds to a calculation in which a single concentration parameter, either 
Co (red), a (green), or (3 (blue), is systematically offset upwards of its fiducial value 
by 10%, while assuming that all of the concentration parameters are known with 
perfect accuracy and precision. The value of the induced systematic error on Wo{ws) 
has been normalized by the statistical constraints on the parameter at the level of 
photo-z calibration specified by the horizontal axis value for Ngpec- We propagate 
systematic errors via Eq. 2.11, as in § 3.3.1. 

The magnitude of the systematic error induced on dark energy parameters mono- 
tonically decreases as Ngpec increases. The physical interpretation of this trend applies 
to nearly all of the results presented in this manuscript, and so we discuss it in de- 
tail in §4. Briefly, as photo-z priors are relaxed the cosmological interpretation of 
the weak lensing signal must rely more heavily on precise knowledge of the matter 
distribution. In the context of the halo model this implies that errors in halo con- 
centrations have more drastic consequences for dark energy parameter inference at 
lower values of Ngpec- 

Dark energy systematics are also less severe when galaxy clustering information is 
included; for each halo parameter, and at every level of photo-z calibration, the solid 
curves are smaller in magnitude than the dashed. This trend is to be expected as its 
analogues have manifested in previous sections: the additional information available 
in galaxy clustering statistics mitigates the consequences for cosmology of errors in 
the matter power spectrum prediction (see Ref. [25] for the analogous benefit of 
mitigating photo-z systematics by including galaxy correlations). Notice that there 
is a proportionally greater mitigation of the systematics at lower values of Ngpec; the 
interpretation of this observation is somewhat subtle. When photo-z priors are weak, 
galaxy clustering information plays a more important role in the self-calibration of 
the photo-z parameters^, a* and z^^^^. As discussed in the preceding paragraph, 
systematic errors induced by incorrect predictions of halo concentrations are more 

^This is an important observation in itself and has been noted elsewhere in the literature (for 
example, Ref. [25]). Most of the so-called "complementarity" of weak lensing and galaxy clustering 
stems from these signals calibrating each other's nuisance parameters. 



-26- 




-2.5 



-3.0 1: 

10^ 



10^ 



10^ 



10^ 



10'' 



10' 



Figure 9: In the top panel we plot the systematic error on wq against Ngpec, the quantity 
encoding the precision with which the photo-z distribution has been calibrated. The same 
quantity is plotted for Wg, in the bottom panel. Dark energy systematics (5(i(;o,w'a) have 
been scaled by (j^p^''^('u;o, if^a), the statistical constraints on the parameter at the level of 
Nspec given by the horizontal axis value of Ngpec- The systematic errors are induced by 10% 
misestimations of our three halo concentration parameters cq, a, and f3. Results for each 
parameter are color-coded according to the legend in the top panel. The concentration 
parameters are defined by c(m, z) = co[m/m=K^o]"(l + z)^. Dashed curves correspond to re- 
sults for an analysis using weak lensing information only, solid curves an analysis employing 
galaxy clustering information together with weak lensing in a joint analysis. 



severe when information about the photo-z distribution is limited. Therefore, because 
including galaxy clustering information at low Nspec has a greater impact on the 
statistical constraints than when Nspec is very large, there is a concomitantly greater 
mitigation of the induced systematics associated with halo concentration errors at 
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lower values of Nspcc- 

In the limit of very large Ngpec, when uncertainty in the photo-z distribution can 
be neglected, a 10% error in the mean halo concentration cq induces a systematic error 
on dark energy parameters that is comparable to or worse than the statistical con- 
straints. This estimation is consistent with previous results (Refs. [28, 82]). For halo 
concentration errors made at other levels, we can estimate the induced systematic 
error on the inferred value of dark energy parameters via simple linear scaling^. For 
example, as can be seen in Fig. 9, in the limit of very large Ngpec a 10% misestimation 
of the parameter cq induces a systematic error 6wq = 1.2a{wQ) for a weak lensing-only 
analysis; thus a 25% error on cq induces a 6wo = [1.2 x (0.25/0.1)] (t(wo) = 3a{wo) 
systematic. For the sake of concreteness we conclude this section with the following 
rough guideline that is based on Fig. 9 evaluated at Ngpec ~ 10^ for a joint analysis: 
in order to guarantee that dark energy systematics induced by halo concentration 
errors are kept at or below the statistical constraints, the parameter cq must be cal- 
ibrated to an accuracy of 5% or better, the parameter a to an accuracy of 12% or 
better, and /3 to better than 25%. 

4. Discussion 

We have studied the significance of matter power spectrum uncertainty for weak lens- 
ing measurements of dark energy. Our results can serve as an updated guideline for 
the calibration requirements on theoretical predictions of Ps{k) and photometric red- 
shift distributions. The photo-z requirements revise those in Ref. [45], who modeled 
nonlinear evolution using the Peacock & Dodds fitting formula, and were thus overly 
pessimistic about photo-z calibration. We also revise the requirements for precision 
in the prediction of Ps{k) outlined in Ref. [33], who assumed perfect knowledge of 
the distribution of photometric redshifts and were thus overly optimistic. 

Both of our models for uncertainty in the theoretical prediction for the matter 
power spectrum have been studied previously. Our second model, in which we allow 
the value of Ps{k) to vary freely about its fiducial values in ten bins of bandpower 
(see §2.1 for a detailed description) is based on the treatment in Ref. [33]. Our 
results are in good agreement with theirs, where applicable. We have generalized 
their results by 1) studying the self-calibration limit ofPs{k) uncertainty, 2) including 
galaxy clustering statistics in the set of observables, and 3) by taking into account 
uncertainty in the distribution of photometric redshifts. Our motivation to treat 
uncertainty in Ps{k) and the photo-z distribution simultaneously comes from results 
presented in the Appendix of Ref. [30], where the authors showed that the photo-z 
calibration requirements vary significantly depending on the assumed fiducial model 

^The stability of the derivatives of our observables with respeet to halo eonccntration parameters 
over a broad range of numerical step sizes ensures the accuracy of this simple linear scaling. 
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of the Ps{k) in the nonhnear regime. This result suggests a nontrivial interplay 
between the photo-z and matter power spectrum calibration demands. 

The cause of this interplay has a simple physical interpretation. In weak lensing, 
there is a degeneracy between the redshift of a galaxy whose image is distorted and 
the typical size of the overdensity responsible for most of its lensing: at fixed angular 
scale, correlations in the image distortions of sources at high redshift are produced 
(on average) by overdensities that are larger in comoving size than those producing 
correlations at low- redshift. The more precisely the photo-z distribution is known, 
the narrower the range of possible wavenumbers contributing to the lensing signal. 
As priors on photo-z parameters are relaxed, the redshifts of the sources are known 
with decreasing precision, and so more information about the power spectrum is 
required in order to compensate. Thus at lower values of Ngpec, dark energy parameter 
inference at a fixed level of statistical uncertainty requires more precise knowledge 
of the matter power spectrum. 

Consequences of this basic physical picture appear throughout this manuscript. 
For example, in Figures 5 and 9, appearing in §3.3.1 and §3.3.3, respectively, the 
degeneracy between source redshift and length scale manifests as dark energy biases 
being more severe for the lower values of Ngpoc- Another example appears in Figure 
4 of §3.2.2, in which we can see that the statistical constraints on wo and degrade 
more rapidly with uncertainty in halo concentrations when Ngpec is small relative to 
larger values of Ngpec. Similarly, the statistical constraints on dark energy discussed 
in §3.2.1 degrade more rapidly as priors are relaxed on 51n(Pi) at lower values of Ngpec. 
The nontrivial relationship between uncertainty in Psik) and in the distribution of 
photometric redshifts clearly illustrates that a detailed and accurate study of the 
calibration requirements on future imaging surveys requires the simultaneous account 
of these contributions to the dark energy error budget that we present here. 

The significance of uncertainty in halo concentrations for the dark energy pro- 
gram has also been studied previously [28, 82]. Again, we have generalized their 
calculations and our conclusions are in good agreement with their results, where 
commensurable. In particular, Ref. [28] studied the prospects for future weak lens- 
ing surveys to self-calibrate halo concentration parameters while simultaneously con- 
straining dark energy parameters. Even with only very modest prior information on 
photo-z parameters, we agree with the conclusion in Ref. [28] that the prospect for 
future imaging surveys to self-calibrate uncertainty in halo concentrations is very 
promising, especially when galaxy correlation statistics are employed in a joint anal- 
ysis: the statistical degradation on wq from self-calibrating the mean, redshift-zero 
halo concentration cq is less than 6% for Ngpec = 5000 (corresponding roughly to 
^o'z/o'z ~ 10~^). The degradation in the constraints is even milder when the photo- 
z distribution is more precisely characterized. However, this result is provisional in 
that it relies on halo concentration being the most significant mode in which the 
power spectrum is uncertain. 



-29- 



Our halo model-based treatment of Ps{k) uncertainty is well- motivated by Ref. [36], 
who carried out a suite of numerical cosmological simulations including hydrodynam- 
ics with a variety of energy feedback mechanisms. One of the salient conclusions of 
Ref. [36] is that the effects of baryonic physics on the matter power spectrum can 
be well-modeled as an enhancement to the mean concentration of dark matter halos. 
Thus we chose our first model of uncertainty in Ps{k) with the intention to study the 
requirements of the dark energy program for precision in our ability to predict the 
effects of baryons on the large-scale distribution of matter. However, recent results 
from the Overwhelmingly Large Simulation (OWLS) project [83] suggest that an 
energy feedback mechanism modeling the effects of AGN is necessary to reproduce 
the characteristics of groups of galaxies [84] . In a recent study based on these results 
[39] , the authors found that a multicomponent halo model with a gas profile that is 
independent from the dark matter profile can accurately model the power spectra 
in the OWLS project. We note, however, that not even N-body simulations have 
achieved the desired precision (< 0.5%) over the full range 0.1 /iMpc^^ ^ k < 5 
/iMpc~^ required of the Psik) calibration [42]. Nonetheless, these results are intrigu- 
ing and suggest that a more complicated model than the one we consider here may 
be necessary to fully encapsulate the baryonic modifications to the matter power 
spectrum. We leave the development and exploration of such a model as a task for 
future work. 

When the assumption that the halo model accurately characterizes all the gross 
features of Psik) is relaxed, self-calibrating the matter power spectrum is very likely 
to be infeasible. To see this, we turn back to our second, more conservative model 
of matter power spectrum uncertainty. As shown in Fig. 2, even in the limit of 
perfectly precise prior knowledge on the photo-z distribution, future surveys are 
unable to self-calibrate the value of Ps{k) to better than 7% on any scale, rendering 
Ps{k) as a dominant component in the error budget and increasing errors on wq 
and Wa_ by a factor of 3 or more. As discussed in §3.3.1, ensuring that dark energy 
systematics induced by misestimations of Ps{k) are kept at or below the level of 
statistical uncertainty of an LSST- or Euclid-like survey, the theoretical prediction 
for the matter power spectrum will need to be accurate to at least 0.5% or better 
on all scales k ^ 5 /iMpc~^. These results reinforce the necessity of an aggressive 
campaign of numerical cosmological simulations if surveys such as LSST or Euclid 
are to achieve their potential as dark energy experiments. 

To contextualize these findings with the current state-of-the-art in numerical 
simulations, we first compare this requirement to the results from the Coyote Universe 
project [41], a suite of nearly 1,000 N-body (gravity only) simulations spanning 38 
fiducial wCDM cosmologies. To date, this is the most ambitious campaign of N-body 
simulations yet performed with the aim to robustly calibrate Ps{k) over the full range 
of scales relevant to weak lensing. The power spectrum emulator based on their 
results has recently been completed [85], and in Ref. [42] the authors demonstrate 
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that results from their simulations can be used to model Ps{k) with sub-percent 
accuracy on all scales /c < 1 /tMpc~^ In addition to the need to expand the range 
of scales over which this level of precision has been attained, the Coyote Universe 
does not account for the nonlinear effects of neutrino mass [86], which have been 
recently established [87] to introduce percent-level changes to the standard Smith et 
al. method of prediction. Moreover, Coyote's 1% precision only applies to wCDM 
cosmologies; extending this level of precision to dynamical dark energy models is an 
active area of current research on in this field, for example, Refs. [44, 88, 89, 90]. 
These complications aside, the Coyote Universe project is, by itself, insufficient to 
completely calibrate the matter power spectrum because N-body simulations neglect 
the effect that baryonic gas has on Ps{k)- Suites of hydrodynamical simulations 
such as the OWLS project [91] discussed above will be essential contributions to the 
calibration program. Continued improvement both in N-body and hydrodynamical 
simulations will clearly be necessary in order to meet the calibration requirements 
we present here. 

4.1 Caveats 

In updating the photo-z precision requirements for future imaging surveys we have 
quantified the uncertainty in the photo-z distribution in terms of Nspec- We reiter- 
ate here an important difference between the meaning of Ngpec in our forecasts and 
elsewhere in the literature. In this work, the quantity Ngpec defines a one-parameter 
family of priors on the photo-z distribution via Eqs. 2.4 & 2.5. In practice the actual 
number of galaxies in the calibration sample will likely need to be larger than Ngpec, 
for example because it will be challenging to obtain a calibration sample that fairly 
represents the color space distribution of the galaxies in the imaging survey. More- 
over, even if such a representative sample is obtained in a particular patch of sky, 
sample variance due to the relatively narrow sky coverage of current and near-future 
calibration samples has a significant impact on the accuracy of the calibration [92] . 
We sought to provide general guidelines for a broad range of future imaging surveys, 
and so we have not attempted to model how these important, survey-specific issues 
affect the calibration requirements. Instead, in our formulation the photo-z preci- 
sion requirements are formally specified in terms of the necessary amount of prior 
knowledge on the photo-z distribution, which in turn is encoded by the parameter 
N 

^ " spec • 

Our modeling of galaxy clustering has several simplifying assumptions that are 
relevant to the calibration of future imaging surveys. In modeling galaxy bias as 
a function of redshift only, we have implicitly assumed perfect knowledge of how 
galaxy bias depends on wavenumber. Uncertainty in the scale-dependence of the 
galaxy bias is degenerate with uncertainty in Ps{k), and so the improvement in the 
constraining power of a survey provided by including galaxy clustering will degrade 
when accounting for uncertainty in the scale-dependence of galaxy bias. However, we 
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have restricted the range of angular scales on which galaxy correlations are exploited 
so that we probe only very large scales (~ SO/i'^Mpc), and so we expect corrections 
accounting for this scale-dependence are small. 

The benefit of galaxy correlations to the photo-z calibration is diluted by lensing 
magnification bias [53, 93], which induces a spurious correlation between sources that 
are well separated in redshift space. This effect thus threatens the ability of galaxy 
cross-correlations to detect and calibrate outliers in the photo-z distribution and will 
need to be accounted for in order to fully realize the potential of galaxy clustering. 
We intend to generalize our results to include these effects in a future paper. 

5. Conclusions 

We have studied the matter power spectrum calibration requirements for future very- 
wide-area weak lensing surveys such as LSST or Euclid. While our findings apply 
to all planned imaging surveys designed to use weak lensing to study dark energy, 
we have phrased our conclusions in terms of these particular surveys because their 
calibration demands are the most stringent. Our results generalize previous find- 
ings by simultaneously accounting for photometric redshift uncertainty, as well as by 
studying the significance of galaxy clustering information. We explored two differ- 
ent models for uncertainty in the nonlinear physics of gravitational collapse, which 
we describe in detail in §2.1. In our first model, we assume that the Halo Model 
accurately predicts the gross shape of Ps{k), but that the internal structure of halos 
is uncertain. In our second, more agnostic model, we allow the value of Psik) to 
vary freely about its fiducial value in ten logarithmically-spaced bins of bandpower 
spanning the range 0.01 /iMpc~^ < A; < 10 /^Mpc"^. We conclude this manuscript by 
providing a brief summary of our primary results. 

1. Future imaging surveys will be unable to self-calibrate the value of Ps{k) to 
better than 7% on any scale. This renders infeasible the possibility of com- 
pletely self- calibrating the theoretical prediction for the matter power spec- 
trum because systematic errors at such levels would induce unacceptably large 
biases in the inferred value of the dark energy equation of state. Moreover, the 
marginalized constraints on Ps{k) are the tightest at scales k ~ 0.2 /iMpc~^, 
nearly an order of magnitude larger in size than where the unmarginalized 
constraints computed in Ref. [33] attain their minimum, emphasizing the ne- 
cessity of a precise calibration of Psi^k) over the full range of wavenumbers 0.1 
hUpc-^ <k<5 hMpc-\ 

2. To ensure that systematics are kept at levels comparable to or below the sta- 
tistical constraints on wq and ifa, Ps{k) must be accurately predicted to a 
precision of 0.5% or better on all scales A; ^ 5 /iMpc"^ in advance of future 
weak lensing observations that will be made by LSST or Euclid. 
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3. The required precision for the Ps{k) prediction as well as the scale to which this 
precision must be attained depend sensitively on iranx, the maximum multipole 
used in the cosmic shear analysis. Figures 7 and 8 together provide a concrete 
guideline that can be used to directly inform the optimal choice for £max that 
balances the need for statistical precision against the threat of matter power 
spectrum systematics. 

4. In keeping with the results in the Appendix of Ref. [30], we find that the 
photo-z calibration requirements are less stringent by a factor of ~ 3 when the 
nonlinear evolution of Ps{k) is modeled with the Smith et al. fitting formula 
relative to Peacock & Dodds, significantly relaxing the demands for photo-z 
precision that appear in Ref. [45]. 

5. Dark energy constraints degrade ~ 40% more slowly with photo-z uncertainty 
when including galaxy correlations in a joint analysis with weak lensing, even 
when the clustering information is restricted to degree-scales and with coarse 
tomographic redshift binning so that Baryonic Acoustic Oscillation features are 
not resolved. 

6. Including galaxy clustering statistics (again, even when Baryon Acoustic Os- 
cillation information is neglected) also significantly relaxes the calibration re- 
quirements and mitigates the severity of systematic errors induced by erroneous 
predictions for Ps{k), especially when prior information on the photo-z distri- 
bution is weak. Dark energy systematics can be reduced by up to 50% by 
including galaxy clustering information; the statistical constraints on wq and 
Wa can degrade 2 — 5 times more rapidly with Ps{k) uncertainty when galaxy 
correlations are neglected. 

7. The redshift-zero, mean halo concentration, cq, must be accurately predicted 
with a precision of 5% or better to keep systematics in dark energy parame- 
ters below the level of statistical constraints. If internal halo structure is the 
dominant mode of Ps{k) uncertainty, then the prospect for self-calibrating cq 
are quite promising, as this would only degrade the dark energy constraints by 
5 - 10%. 

8. The matter power spectrum calibration requirements are more stringent when 
the distribution of photometric redshifts is known with less precision. This 
effect is due to a degeneracy between source redshift and lens size, and is the 
chief motivation for a simultaneous account of these sources of uncertainty. 
We find that the constraints on wq and degrade 2 — 3 times more rapidly 
with Ps{k) uncertainty for spectroscopic calibration samples with the statistical 
equivalent of Nspec ~ 5000 relative to Ngpcc ~ 10^. 
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9. The requirements for the precision with which Ps{k) need be predicted are, in 
general, less stringent for DES than for LSST by a factor of a few. To ensure 
that matter power spectrum systematics do not contribute significantly to the 
dark energy error budget for DES, we find that if correlations up to a maximum 
multipole of ^max = 3000 are used in the lensing analysis then Ps{k) will need to 
be calibrated to an accuracy of 2% or better on scales k < 5 /iMpc^^. The DES 
requirements scale with ^^ax in a similar fashion to the scaling of the LSST 
requirements summarized in Fig. 8. The difference between the requirements 
is driven by the relative depth of these two surveys: the constraining power 
on dark energy provided by cosmic shear measurements derives chiefly from 
small-scale correlations {k > l/iMpc"^) where shot noise is most significant. 
The shallower depth of DES (^med = 0.7, A^a = 15 gal/arcmin^) results in 
these modes being less informative, and so DES suffers less from uncertainty 
in small-scale information. 
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